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Within the weak-field, post-Newtonian approximation of the metric theories of gravity, we deter- 
mine the one-way time transfer up to the order 1/c'', the unperturbed term being of order 1/c, and 
the frequency shift up to the order 1/c''. We adapt the method of the world-function developed by 
Synge to the Nordvedt-Will PPN formalism. We get an integral expression for the world-function 
up to the order l/c^ and we apply this result to the field of an isolated, axisymmetric rotating body. 
We give a new procedure enabling to calculate the influence of the mass and spin multipole moments 
of the body on the time transfer and the frequency shift up to the order 1/c'*. We obtain explicit 
formulas for the contributions of the mass, of the quadrupole moment and of the intrinsic angular 
momentum. In the case where the only PPN parameters different from zero are /3 and 7, we deduce 
from these results the complete expression of the frequency shift up to the order 1/c'*. We briefly 
discuss the influence of the quadrupole moment and of the rotation of the Earth on the frequency 
shifts in the ACES mission. 



I. INTRODUCTION 



Owing to recent progress in absolute frequency measurements of some optical transitions with a femtosecond 
laser, it seems possible to achieve in a near future atomic clocks having a time-keeping accuracy of the order of 
10"^* in fractional frequency [Q. Since a reduced gravity would significantly increase clock performances, it is 
envisaged to install such clocks on board artificial satellites and to compare them with terrestrial clocks by exchange 
of electromagnetic signals. Already, a spatial experiment like ESA's Atomic Clock Ensemble in Space (ACES) mission 
^ is planned for 2006, the purpose being to obtain an accuracy of order 10~^^ in fractional frequency. 
At a level of uncertainty about 10~^^, a fully relativistic treatment of time/frequency transfers must be performed 



up to the order 1/c'* | 
transfer, the main re 



23( 1 . As far as we know, the corresponding calculations have not been carried out. For the time 
ativistic correction of order 1/c'^ is the well-known Shapiro time delay Q]. Other corrections 
due to the quadrupole moment and to the intrinsic angular momentum have been studied by several authors 
Gravitational corrections of order 1/c^ in the frequency transfers were theoretically determined and experimentally 
checked a long time ago These corrections are now commonplace in the Global Positioning System. The relativistic 
theory of the frequency transfers have been recently extended up to the terms of order 1/c'^ 0, justifying the results 
previously given in § without any detail. However, it must be pointed out that, in 0, some terms of order l/c"^ 
due to the quadrupole moment J2 of the Earth are bounded without any explicit calculation. Furthermore, the 
time/frequency transfers have been calculated within the limited framework of general relativity, which prevents from 
discussing new tests of gravitational theories. 

The present paper is a first step towards a general theory of the time and frequency transfers including all the terms 
of order 1/c^, within the post-Newtonian approximation of any metric theory of gravity. Using the Nordtvedt-Will 
PPN formalism we bring the complete determination of these effects in the field of an isolated, axisymmetric 
rotating body, the gravitational field being assumed stationary. Of course, modelling a mission in the vicinity of the 
Earth at a level of accuracy about 10~^* will require to add the effects due to the tidal gravitational field induced by 
the Sun and the Moon. 

We assume that the photons ensuring the transfers follow null geodesies. The problems that we have to tackle come 
down to the following ones, relative to a couple of points xa = {ctA,XA) and xb = {ctB^XB) connected by a null 
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geodesic: 

i) to calculate the (coordinate) time transfer tg — as a function of (a;^, Xg); 

a) to determine the vectors tangent to the null geodesic at xa and xb- Solving this second problem is indeed 
indispensable to calculate the frequency shift between xa and xb- 

The method generally employed to study the questions related to the p rop agation of light in a gravitational field 
is based on the solution of the null geodesic equations (see, e.g., H, |ll|^ |l^, [ij for investigations in the linearized, 
weak- field limit of general relativity). However, the theory of the world- function developed by Synge jl^ presents the 
great advantage to spare the trouble of integrating the geodesic equations. Once the world-function is determined, it 
is possible to solve straightforwardly the two above-mentioned problems. This method is particularly elegant for the 
stationary, axisymmetric field and we apply it in the present paper. We find a new procedure enabling to determine 
the influence of the mass and spin multipole moments of the body. Explicit calculations are given for the contributions 
of the mass, of the quadrupole moment and of the intrinsic angular momentum of the rotating body. 

The paper is organized as follows. In Sec. II, the relevant properties of the world-function ^{xa, xb) are recalled and 
the general expression of this function in the post-Newtonian limit of any metric theory is given. The corresponding 
expression of the time transfer tB — tA is derived up to the order 1/c^. In Sec. Ill, we determine the expression 
of Q,{xa,xb) and of ts — tA within the ten-parameter PN formalism of Nordtvedt and Will. Then, in Sec. IV, we 
focus our attention on the case of an isolated, axisymmetric rotating body. We show that it is possible to determine 
the contributions of the mass and spin multipole moments by straightforward differentiations of a single function. 
Retaining only the terms due to the mass Af , to the quadrupole moment J2 and to the intrinsic angular momentum 
S of the rotating body, we obtain explicit expressions for the time transfer up to the order and for the tangent 
vectors at xa and xb up to the order 1/c^. In Sec. V, the frequency shift is developped up to the order in the 
case where /3 and 7 are the only non vanishing PPN parameters. We find detailed expressions for the contributions 
of J2 and S and we discuss the possible influence of these terms in the ACES mission. We give our conclusions in 
Sec. VI. 

In this paper, G is the Newtonian gravitational constant and c is the speed of light in a vacuum. The Lorentzian 

metric of space-time is denoted by g. The signature adopted for 5 is (-1 ). We suppose that the space-time is 

covered by one global coordinate system (x'') = (x", a;), where x^ = ct, t being a time coordinate, and x = (x*), the x* 
being quasi Cartesian coordinates. We assume that the curves of equations x' = const, are timelike, which means that 
goQ > anywhere. We employ the vector notation a in order to denote either (a^ , o^, a'^) = (a*) or (ai, 02, 03) = (a^). 
Considering two such quantities a and b with for instance a = (a*), we use a • 6 to denote a*6* if b = (6*) or a^bi if 
b = (bi) (the Einstein convention on the repeated indices is used). The quantity | a \ stands for the ordinary Euclidean 
norm of a. 



II. THE WORLD-FUNCTION AND ITS POST-NEWTONIAN LIMIT 
A. Definition and fundamental properties 

Consider two points xa and xs in a space-time with a given metric g^i, and assume that xa and x b sue connected 
by a unique geodesic path F. Throughout this paper, A denotes the unique afhne parameter along P which fulfills the 
boundary conditions = and As = 1. The so-called world-function of space-time is the two-point function 
il{xA,XB) defined by 

n{xA,XB) = -J^ 3,Ax''{>^))^-^^d\^ (1) 

the integral being taken along P. It is easily seen that VI{xaiXb) — e\sAB^ 1"^, where sab is the geodesic distance 
between xa and xb and e = 1,0,-1 for timelike, null and spacelike geodesies, respectively. It results from definition 
(|l|) that the world-function f2(xA,a;s) is unchanged if we perform any admissible coordinate transformation. 

The utility of the world- function for our purpose comes from the following properties |l4| . 

I) Two points Xa and xb are linked by a light ray if and only if the condition 

r!(x^,XB)=0 (2) 

is fulfilled. Thus, f2(xA, x) = is the equation of the light cone C(x a) at xa- This fundamental property shows that if 
r2(xA, a^s) is known, it is possible to determine the travel time tB — tA of a photon connecting two points xa and xb 
as a function of tA, xa and xb- It must be pointed out, however, that solving the equation ^l{ctA^XA,ctB,XB) = 
for tB yields two distinct solutions t~^ and t^ since the timelike curve x* = x^ cuts the light cone C{xa) at two points 
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Xq and Xq, x~^ being in the future of x~^. In the present paper, we always regard xa as the point of emission of the 
photon and as the point of reception, and we are concerned only with the determination of tg — i^i as a function 
of tA, Xa and xb- We put 

t^~tA^T{tA,XA,XB), (3) 

and we call T{tA, xa, xb) the (coordinate) time transfer function. Of course, it is also possible to introduce an other 
time transfer function giving — as a function of the instant of reception and of a; a, xb^ but we do not use it 
here. 

ii) The vectors (dx°' /dX)A and {dx"/dX)B tangent to the geodesic F respectively at xa and xb are given by 



lap 



dx/^N an / dxP\ _ an 



dX J dx"^ \ dX J g dx'^ 



(4) 



As a consequence, if n{xA, xb) is explicitly known, the determination of these vectors does not require the integration 
of the differential equations of the geodesic. Let us note that it can be proved that the tangent vectors (^) are null 
when (H) holds. 

Consider now a stationary space-time. In this case, we use exclusively coordinates {x'^) such that the metric does 
not depend on x . Then, the world-function is £i function of — x^^ and Xb, and ^ reduces to a relation of 
the form 

tB-tA^T{xA,XB). (5) 

The time transfer function T{xa,xb) plays a central role in the present paper because a comparaison between (|^) 
and (^) immediately shows that the vectors {1'^)a and {1'^)b defined by their covariant components 

d 

{Io)a = 1, {k)A = c---^T{xa,xb) , (6) 
ox^ 

d 

{h)B^l, {h)B = ~c-—T{xa,xb) , (7) 

OXb 

are tangent to the ray at xa and xb, respectively. It must be pointed out that these tangent vectors correspond to 
an affine parameter such that Iq = I along the ray (note that such a parameter does not coincide with A). Generally, 
extracting the time transfer formula (0) from (|^), next using (^-(0) will be more straightforward than deriving the 
vectors tangent at xa and xb from (|j)J next imposing the constraint (||). Wc shall use (||) and (|^) in Sec. IV. 

To conclude, let us emphasize that the method of the world-function works as long as Q{xa, xb) is a well-defined, 
single- valued function of xa and xb- This condition is satisfied in any region of space-time in which any points xa 
and Xb are connected by one and only one geodesic, a feature which excludes the existence of conjugate points. This 
requirement is certainly fulfilled in experiments performed in the solar system and more generally for observations of 
stars belonging to our Galaxy. 

B. General expression of the world- function in the post-Newtonian limit 

To begin, let us assume that the metric may be written as 

gfj^L' = Vf^i^ + h^i, (8) 

throughout space-time, with 77^^ = diag(l, —1, ^1, —1). Let F^q) be the straight line defined by the parametric 
equations = x"q^{X), with 

xf,){X) = {x%-x'X)X + x'X, 0<A<1. (9) 
With this definition, the parametric equations of the geodesic F connecting xa and xb may be written in the form 

x"(A)=a;^o)(A)-fX"(A), 0<A<1, (10) 
where the quantities X"(A) satisfy the boundary conditions 

X"(0) = 0, X"(1) = 0. (11) 
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Inserting and dx^{\)/dX = — + dX^(X)/dX in (|^), then developing and noting that 

„ . dX" 



dX 



-dX = 



by virtue of (Jl 1|) , we find the rigorous formula 

1 

n{xA,XB) = n^''HxA,XB) + M~x'X)ix''B~-x'X) h^^ix'^ {X))dX 



2 7o 



dX^ dX" 



(12) 



(13) 



dX dX 



+ 2(x^-x:^)V(^"(A)) 



dX"" 



dX 



dX, 



where the integrals are taken over F and n^^\xA, a^s) is the world-function in Minkowski space-time 

fl^°'>{xA,XB) = ^Vf^i'ixB - ^Di^B - X^a) . 



(14) 



Henceforth, we shall only consider weak gravitational fields generated by self-gravitating extended bodies within 
the slow-motion, post-Newtonian approximation. So, we assume that the potentials h^i, may be expanded as follows 



"OO — ""oo "I ^'*n 



^4 "00 



0(6). 



0(4). 



From these expansions and from the Euler-Lagrange equations satisfied by any geodesic F, namely 



dx^^ 



1 



dx^ dx" 



dX [^"^ dX J " 2^'''''"' dX dX ■ 



(15) 



(16) 



it results that X''(A) = 0(2) and that dxf'/dX = x'^ -x'^ + 0{2). As a consequence, h^^ix^'iX)) = /i^^(x^p)(A))-|-0(4) 
and the third and fourth terms in the r.h.s. of Eq. ( p^ ) are of order 1/c^. These features result in an expression for 
fl{xATXB) as follows 



where 



n{xA,XB) = n^°\xA,XB) + f7(^^)(x^,XB) + 0(4) , 



4'o^(xfo)(A))dA 



(17) 





Jo 




- XA)iXg - 


A(^B- 





(18) 



the integral being taken over the line F(o) defined by (^). 

The formulas ( p7| ) and (|l8|) yield the general expression of the world- function up to the order within the 
framework of the 1 PN approximation. We shall see in the next subsection that this approximation is sufficient to 
determine the time transfer function T{tA, xa, xb) up to the order It is worthy of note that the method used 

above would as well lead to the expression of the world-function in the linearized weak-field limit previously found by 
Synge 

We s: 



lall put henceforth Rab = xb — xa and Rab =| Rab \- Defining the quantities = (x^ — x^)/Rab, (|lj) 



and (|1S) might be easily rewritten with these notations. 

C. Time transfer at the order l/c^ 

Suppose that xs is the point of reception of a photon emitted at xa- Taking ( p7[ ) into account, (^ may be written 
in the form 

r!(°)(x^,XB) -|-f](^^)(xA,XB) = 0(4) , 
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which imphes the relation 



c cHab 
Using iteratively this relation, we find for the time transfer function 

T-/, X 1 D n'^^^^ctA,XA,ctA + Rab,xb) , 

T{tA,XA,XB) = -Rab ^ + 0(5) . 

c cRab 



(19) 



(20) 



This formula shows that the time transfer T(t^, a;^, a;^) can be explicitly calculated up to the order 1/c^ when 
,xb) is known. This fundamental result will be exploited in the following sections. 
The quantity ^^^'^\ctA,XA,ctA + Rab,xb) in ( pO| ) may be written in an integral form using (p^, in which Rab 
and RabX + ctA are substituted for — x'a and for a;|'Q-|(A), respectively. Hence 



T{tA,XA,XB) - -Rab f' ko^(^"(A)) + h[f{z^{X))N^N^ + -^f (z"(A))iV^ 

the integrals being taken over the line defined by the parametric equations x" = z°'{X), where 

z"(A) = i?^sA + ciA, z'{X) RabN'X + x'a, 0<A<1. 



0(5) 



(21) 



(22) 



It must be noted that the line defined by (E2|) is the null geodesic of Minkowski metric from xa, the direction cosines 



of which are N"^ 



x\)/R 



ab- 



III. WORLD-FUNCTION AND TIME TRANSFER WITHIN THE NORDTVEDT-WILL PPN 

FORMALISM 



A. Metric in the 1 PN approximation 

In this section, we use the Nordvedt-Will post-Newtonian formalism involving ten parameters /3, 7, ^, ai, . . ., ^4 
We introduce slightly modified notations in order to be closed of the formalism recently proposed by Klioner and 
Soffel as an extension of the post-Newtonian framework elaborated by Damour, Soffel and Xu |17 for general 
relativity. In particular, we denote by the velocity of the center of mass O relative to the universe rest frame [ p4[ . 

Although our method is not confined to any particular assumption on the matter, we suppose here that each source 
of the field is described by the energy-momentum tensor of a perfect fluid 



1 + 4 (n + - 

c \ P 



P9 



where p is the rest mass density, 11 is the specific energy density (ratio of internal energy density to rest mass density), 
p is the pressure and is the unit 4-velocity of the fluid. In this section and in the following one, v is the coordinate 
velocity dx/dt of an element of the fluid. We introduce the conserved mass density p* given by 



1 



1 /I 



P = = P 

where g = det{g^^) and U is the Newtonian- like potential 

U{x\x)=G ( P^^^^d'x' 



+ 37C/ + 0(4) 



(23) 



(24) 



In order to obtain a more simple form than the usual one for the potentials ft-oi, we suppose that the chosen (x^) 
are related to a standard post-Newtonian gauge (x^) by the transformation 



1 



= + ^ [(1 -f 2e + ^2 - Ci)dtX - 'io^2Vr ■ Vx] , x' 
where x is the superpotential defined by 

x{x\x)^]^G J p*ix",x') Ix-x'ld'^x'. 



(25) 



(26) 
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Moreover, we define p by 
P = P 



1 + ^(27 + 1 - 2e + as + Ci)^ + (1 - 2/5 + e + C2) J + (1 + Cs)^ + (87 - 2^ + 3C4) ^ 



p*c?- 



-{ai - a3)-7 - - 2a3)—: ^ + 0(4) 



Then, the post-Newtonian potentials read 



where 



2/5 2 2^ 



2a2 



-^w + ^u.^ + -^cj^w + ^(Ci - 200. - ^z;>^a,,x + 0(6), 



h = {Hq^} = — 



1 \ 1 
7 + 1 + -ai j w + -aiw V, 



0(5), 



27 



w6^j + 0(4) 



2 



/9*(a;0,a;')p*(a;°,a;")(a;-a;') f x' - x" x - x" 



\x~x'\'^ \\x — x" \ \X' — X" 



(fx'(fx'\ 



I 13 



w{x°,x) j 



\X — X' ^ 



(27) 

(28) 
(29) 
(30) 

(31) 
(32) 
(33) 
(34) 



B. Determination of the world-function and of the time transfer 



For the post-Newtonian metric given by (g8|)-(^4|), it follows from ( |18[ ) that Vl{xatXb) niay be written up to the 
order l/c"^ in the form given by Eq. (17) with 



{XA ,Xb) = n^r^ [XA ,Xb) + ^) (XA ,Xb) + f^i:^) (XA , xb) 



where 



n<-S^\xA,XB) 



«;(a;fo)(A))dA, 

nij^""^ {xA,XB) = ^(^^ + l + ^ai^ {x% - x"a)Rab ■ w{xf„) {X))dX , 
ili^^^ {xA, xb) = {x% - x\){Rab ■ Vr) w{xf,^ {X))dX , 



the integrals being calculated along the line defined by 



The corresponding time transfer function is easily obtained by using (EG) or (21). We get 



T{tA,XA,XB) = -Rab + \{i+1)Rab I w{z'^{X))dX 
c c Jq 



2 „ 

-^Rab 



(7+1 + iai) / w{z"{Xj)dX + ]aiVr I w{z"{X))dX 



0(5), 



(35) 

(36) 
(37) 
(38) 



(39) 



the integral being evaluated along the curve defined by (|2£ 

Let us emphasize that, since w = U + 0(2), w may be replaced by the Newtonian-like potential U in expressions 
(§)-(§). 
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C. Case of stationary sources 



In what follows, we suppose that the gravitational field is generated by a unique stationary source. Then, dtX ~ 
and the potentials w and w do not depend on time. In this case, the integration involved in (|36|)-(|3^) can be performed 
by a method due to Buchdahl [^. Introducing the auxiliary variables ua = — x' and ys Xb~x' , and replacing 
in ^ the parameter A by m = A — 1/2, a straightforward calculation yields 

w{x^a){X))d\ = G [ p{x')F{x',XA,XB)d^x', (40) 



^ w{x^o)W)dX^G J p*{x')v{x')F{x\xA,XB)d^x' , 



(41) 



where the kernel function F{x' , xa, xb) has the expression 

"1/2 du 



F{x,Xa,Xb)= r —r, \ ^ 

J -1/2 I (.Vb - Va)u + ^{vb + Va) 



Noting that Ub — Ua — Rab, which implies that Iub — yA\= Rab, we find 

r^t N 1 1 (\X-Xa\ + \X-Xb\+Rab\ 

F{x,xa,xb) = - — In I ^ — . 42) 

Rab \\x-xa\ + \x-xb\-Rab J 

Inserting (^), (^) and ( p^ ) in (^)-(^) and in (|3^) will enable one to obtain quite elegant expressions for 
^'^p^){^xa,xb) and for T{xatXb), respectively. 



IV. ISOLATED, AXISYMMETRIC ROTATING BODY 



Henceforth, we suppose that the light is propagating in the gravitational field of an isolated, axisymmetric rotating 
body. The gravitational field is assumed to be stationary. The main purpose of this section is to determine the 
influence of the mass and spin multipole moments of the rotating body on the coordinate time transfer and on the 
direction of light rays. From these results, it will be possible to obtain a relativistic modelling of the one-way time 
transfers and frequency shifts up to the order 1/c'* in a geocentric non rotating frame. 

Since we treat the case of a body located very far from the other bodies of the universe, the global coordinate 
system (x^) used until now can be considered as a local (i.e. geocentric) one. So, in agreement with the UAI/UGG 
Resolution Bl (2000) we shall henceforth denote by W and W the quantities w and w respectively defined by 
( ^ ) and (jsj) and we shall denote by G^i/ the components of the metric. However, we shall continue here with using 
lower case letters for the geocentric coordinates in order to avoid too heavy notations. 

The center of mass O of the rotating body being taken as the origin of the quasi Cartesian coordinates (a;), we 
choose the axis of symmetry as the x'^-axis. We assume that the body is slowly rotating about Ox^ with a constant 
angular velocity a;, so that 

v{x) ^ LJ X x . (43) 

In what follows, we put r =| x |, =| xa \ and vb =\ xa \- We call 9 the angle between x and Ox^. We consider 
only the case where all points of the segment joining xa and xb are outside the body. We denote by the radius 
of the smallest sphere centered on O and containing the body (for celestial bodies, re is the equatorial radius). In 
this section, we assume the convergence of the multipole expansions formally derived below at any point outside the 
body, even if r < r^. 



A. Multipole developments of W and W 



According to (|3l|), (|3^ ) and (^3|), the gravitational potentials W and W obey the equations 

V^IF = -inGp, V^W = -47rGp*a; x x . 



(44) 
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It follows from ( p4[ ) that the potential W is a harmonic function outside the rotating body. As a consequence, W 
may be expanded in a multipole series of the form 



W{x) 



GM 



OO 

X — > / T \ ^ 



n=2 



(45) 



In this development, the Pn are the Legendre polynomials and the quantities M, J2 . . . J„ . . . correspond to the 
generalized Blanchet-Damour mass multipole moments in general relativity [ p^ . 
In fact, taking into account the identity 



9z" \ r 



(-l)"n! 



Pniz/r), z^x^ 



it will be much more convenient for the computation of integral ( ^0| ) to use the following expansion in a series of 
derivatives of l/r 



^ 00 



(-11: 7, 
r — ' n\ oz 

n—2 



W{x)= GM 

According to (^) , the mass density p can be developped in the multipole series 

p{x) = M 



(46) 



(47) 



5^'^\x) being the Dirac distribution supported by the origin O. 
Now, substituting (Esl) into (p4) yields for the vector potential W 



W{x) ^G j 



p*{x')ui X x' 



\x — X' 



It is possible to show that this vector may be written as 



W 



X VV, 



(48) 



(49) 



where V is an axisymmetric function satisfying the Laplace equation V^V = outside the body. Consequently, we 
can expand V in a series of the form 



V{x) = ^ 
r 



00 

P„(cos( 



(50) 



where / and the Kn are constants. Inserting (|5C|) into (^9|) and using the identity 

{n + l)P„(z/r) + (z/r)P;,(z/r) = PL+ii^/r) , 
we find for W an expansion as follows 

GIu) X X 



Wix) 



2r3 



(51) 



which coincides with a result previously obtained by one of us po| . The coincidence shows that / is the moment of 
inertia of the body about the z-axis. Thus, the quantity S = lu) is the intrinsic angular momentum of the rotating 
body. The coefhcients Kn are completely determined by the density distribution p* and by the shape of the body 
pO|, Expansion (M) may also be written as 



W{x) 



1 



GS' X V 



^ 00 



, ^ I ear) 

n—l 



Consequently, the density of mass current can be developped in the multipole series 



p*{x){uj X x) 



1 



-S' X V 



Qn 



dz"^ 



(52) 



(53) 



a property which will be exploited in the following subsection. 



9 



B. Multipole structure of the world-function 



The function fl^^^'^ {xa, xb) is determined by (p5|)-(|38|) where w and w are respectively replaced by W and W. 
The integrals involved in the r.h.s. of (^5|)-(^8|) are given by @ and (|l]). Substituting Q into (|o|) and using the 
properties of the Dirac distribution, we obtain 



T^(a;(o)(A)) dX^GM 
Similarly, substituting (|3|) into we get 



oo 



gn 

9i" 



F(a;,a;A,a;B) 



(54) 



x=0 



W{x^o)W)dX = --GSxV 



°° 1 an 

1- V-if r" — 



F{x,xa,xb] 



(55) 



x=0 



These formulas show that the multipole expansion of fj'^^' {xa, xb) can be thoroughly calculated by straightforward 
differentiations of the kernel function F{x,xa,xb) given by (^2|). They constitute the essential result of the present 
paper, from which it would be possible to deduce the multipole expansions giving the time transfer and the frequency 
shift between A and B up to the order 1/c'*. 

In order to obtain explicit formulas, we shall only retain the contributions due to M, J2 and S in the expansion 
yielding ^[^^"^ and fi^^"*. Then, denoting the unit vector along the z-axis by fe and noting that S = Sk, we get for 

^wixA,XB) 



^w^\xa,xb) 



CM {x% -x^f^ jR\b ,__(rA + rB+ Rab 



2GM 
GM 



Rai 



■In 



\rA + rB - Rai 



[{rA + rBr~R ^ 



J2r^ '-^^ 



ABl 



TA 



tb 



(56) 



[rA+TB Y - Rab 



(k X xa^ (fc X xb)^ 



and for Q^^^\xa, xb) 



(PN) 
W 



(xa,xb) 



7 + 1 + -ai 



2GS 



(A 



^ rA+TB k-jxA^-XB) 
TATB (ta + Tb)"^ - R\b 



(57) 



Finally, owing to the limit | ai |< 0.02 furnished in wc shall henceforth neglect all the multipole contributions in 
^^■uI'^\xatXb)- Thus, we get 



^Z'\xa,xb) = a,f^ix% - .^\)^^ln (!iL±l£±|i£ 
2c'^ Rab \rA + rB- Rab 



(58) 



In this section and in the following one, the symbol + • • • stands for the contributions of higher multipole moments 
which are neglected. For the sake of brevity, when + • • • is used, we systematically omit to mention the symbol 0{n) 
which stands for the neglected post-Newtonian terms. 



C. Time transfer function up to the order 



Let us substitute Rab for x^ — x^ into (|5^)-(^8|) and insert the corresponding expression of fi^^^' into (^ 
get an expression for the time transfer function as follows 



We 



T{xa,xb) = -Rab +Tm(xa,xb) + Tj^{xa,xb) + Ts{xa,xb) + %^{xa,xb) 
c 



(59) 



where 



n- ( \ t ( ta + tb+Rab 

Tm{xa,xb) = (7 + 1)^^ In ■ — ,, 

\rA + rB - Rab J 



(60) 
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Tj^{xa,Xb) = 
Ts{xa,Xb) = 



- 7 + 1 — 

/^0 



J2rlRAB 



[rA+rsY-R'^A 



AB 



2{rA+rB) 
{rA + rsY-R'^ 



AB 



rA 



rB 



{k X xaY X xb) 



. , ^ 1 \2GSta+ rB 
— 7 + 1 + -ai 



k ■ (xA X Xb) 



^'^ tatb {rA+rBY 

GMRAB-Vr, fvA+rB + RAl 



r?2 ■ 

^AB 



2c4 



Rab 



\rA +rB - Rai 



(61) 
(62) 
(63) 



The time transfer is thus expHcitly determined up to the order l/c**. The term of order 1/c^ given by (|60| ) is 
the weh-known Shapiro time delay. Equations ( |6l| ) and ( |6^ ) extend results previously found for 7 = 1 and ai = 
1^. However, our derivation is more straightforward and yields formulas which are more convenient to calculate the 
frequency shifts. As a final remark, it is worthy of note that T/v/ and Tj^ are symmetric in {xa,xb), while Ts and 
are antisymmetric in {xa,xb)- 



D. Directions of light rays at xa and xb up to the order 1/c 



In order to determine the vectors tangent to the ray path at xa and xs, we use Eqs. (0) and (|^) where T is replaced 
by the expression given by (p9|). For the sake of brevity, we put henceforth I a — {{li)A} and Ib = {{li)B}- We find 



lAixA,XB) = -Nab +Im{xa,Xb) +Ij^{xa,Xb) +lsixA,XB) +lv,ixA,XB) H 

Ib{xa,xb) = -Nab - Im{xb,xa) - Ij^{xb,xa) + Is{xb,xa) + 1„Xxb,xa) H 

where Zm, Ij2i and stand for the contributions of 7j\,/, Tj^, Ts and 7^^, respectively. Putting 



we get from (|60| 



From (|T|), we get 



Xa Xb - r - Xa 

nA = — , nB = — , Nab = — 5 

rA rB Rab 



, , . , ^ ^^2GM {rA+rB)NAB + RabUa 
Im(xa,xb) = -(7 + 1) — 2 



(rA + reY-RlB 



IjJxa^xb) = (7+ 1) 2 



rA + rB 



[{rA + rBY-RU' 



xiNAi 



2{k ■ riA + k ■ ub 



2 {rA+rBY + SR\B 
{rA + rsY-R'AB 



rB 



rA + rB 



+2nA- 



l-jk-UA Y , l-ik-nBY\ {rA+rBY + R 
rA 
Rab 
rA + rB 



'-AB 



(rA + rBY - Rab 



1 



1 — 3(fe • tia) 



{3rA + rB){rA + rB) - R 



AB 



(64) 
(65) 



(66) 



(67) 



+ {rA + rB) 
Rab 



2{k ■ nA){k ■ ub) 1 — (fc • tib) 



rA 



rB 



-4fc- 



rA 



(fc . r^^)i^^A+rB){rA+rB) -R'AB ^ . 



2rA{rA + ?'b) 

From (|6^) and (63), we derive the other contributions which are of order l/c" 



IsixA, Xb) = ( 7 + 1 + 



2GS 



rA + rB 



c3 rA[{rA + rBY-R\B] 



k X ub 
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2rArB_k ■ (n^ x n^) 
GM 



{^rA+rB){rA +rB) 



R2 
^AB 



ai- 



Vr - {vr ■ Nab) Nab , / rA 
In 



'2rA{rA + tb) 

rB + Rab 



riA + riB 



2Rai 



+ {vr ■ Nab) 



ta + tb- Rab 

{ta + rB)NAB + RABriA 



{rA + rB)^-R'AB 



(68) 



(69) 



We note that the mass and the quadrupole moment yield contributions of order 1/c^, while the intrinsic angular 
momentum and the velocity relative to the universe rest frame yield contributions of order 1/c^. 

E. Sagnac terms in the time transfer function 

In an experiment like ACES, recording the time of emission tA will be more practical than recording the time of 
reception Ib- So, it will be very convenient to form the expression of the time transfer T{xa,xb) from a;^(i^) to 
XB{t^) in terms of the position of the receiver B at the time of emission tA- For any quantity QB{t) defined along 

the world line of the station B, let us put Qb — Q{tA)- Thus we may write XB{tA), ^B{tA), VB{tA), vb =\vb |, etc. 

Now, let us introduce the instantaneous coordinate distance Dab — xb — xa and its norm Dab- Since we want 
to know Ib — tA up to the order \/c^^ we can use the Taylor expansion of Rab 



R 



AB 



D 



AB 



it 



B - tA)VB 



1, 

-,itB 



tA) aB 



-^{tB - tA)"^ bB 



where as is the acceleration of B and = daB /dt. Using iteratively this expansion together with (pOj), we get 



T{xA, xb) = T{xA, xb) + -^Dab ■ VB 



1 
2? 



Dab 



{Dab ■ Vb) 

^AB 



v 



+ 



1 

1 D 



{Dab ■ Vb) (w| + Dab ■ as) 



:^AB 



B ^ Dab ■ o-B 
1 



Vb ■ aB 



-D 



AB 



JB 



(70) 



AB 



c Dab 



Vb [Tm{xa,xb) + Tj2{xa,xb)] 



\dabVb ■ [Im{xb,xa) + l.h {xb, xa)] 



where T{xa, xb) is obtained by substituting xb^tb and Dab respectively for xb, tb and Rab into the time transfer 
function defined by (^9|)-(|63|). This expression extends the previous formula to the next order The second, 

the third and the fourth terms in ( [70| ) represent pure Sagnac terms of order 1/c^, l/c^ and 1/c^, respectively. The 
fifth and the sixth terms arc contributions of the gravitational field mixed with the coordinate velocity of the receiving 
station. Since these last two terms are of order 1/c'*, they may be calculated for the arguments {xb,xa) (note the 
order of the arguments in Im and Zjj). 



V. FREQUENCY SHIFT IN THE FIELD OF AN AXISYMMETRIC ROTATING BODY 

A. General formulas 



Consider a clock Oa on A and a clock Ob on B delivering respectively the proper frequencies /a and /b and 
suppose that Oa is sending photons to Cb. The one-way frequency transfer from Oa and Ob is characterized by 
the ratio /a/Zb which may be written as Ja/Zb = {fA/vA){vA/vB){vB/fB) where va is the proper frequency of the 
photon as measured on A at the instant of emission and i/b is the proper frequency of the same photon as measured 
on B at the instant of receipt. The ratios /a/i^a and ]bIvb are obtained by local measurements performed on A 
and B, respectively Q. So, in the present paper, we are concerned only with the theoretical determination of va/vb- 
This ratio is given by the well-known relation 

^ = (71) 

>^B Wg{l^)B 
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where = {dx^ /ds)A and = {dx^ /ds)B are respectively the unit 4-velocity of the clock Oa and of the clock 
and {1^)a and {1^)b are the null tangent vectors at the point of emission x^ and at the point of reception XB^ 
respectively. 

Let us denote by va = {dx/dt)A and vb = {dx/dt)B the coordinate velocities of the clocks on A and B, respectively. 
Since the gravitational field is assumed to be stationary, the formula (|7^) giving the frequency shift between xa and 
xb may be written as 



i^A 

VB 



<IA 



qA = l + -Ia - va, qb = I + -Ib ■ vb , (72) 

Qb c c 

where I a and Ib are the quantities respectively defined by (^) and (0). 

It is possible to calculate the ratio qa/qb up to the order 1/c* from our results in Sec. IV since I a and Ib are given 
up to the order l/c^^ respectively by (Q) and (|65|). Denoting by l'-"^ /c" the 0{n) terms in I, qa/qb may be expanded 

as 



QB 



1 



1 Nab ■ {vA 

c 



Vb) 



1 



l-AT, 



Vb 



AB 



(2) 



Va 



7(2) 

Ib ■ '"B 



1 r 



;(3) 7(3) 

Ia ■ ""A - Ib ' '"b 



AB 



{if ■ Vb) {va - 2vb) + {if ■ va) Vb] + 0(5) 



(73) 



In order to be consistent with this expansion, we have to perform the calculation of u\/vPg at the same level of 
approximation. For a clock delivering a proper time r, \/vP is the ratio of the proper time dr to the coordinate time 
dt. To reach the suitable accuracy, it is therefore necessary to take into account the terms of order 1/c* in goo- For 
the sake of simplicity, we shall henceforth confine ourselves to the fully conservative metric theories of gravity without 
preferred location effects, in which all the PPN parameters vanish except (3 and 7. Since the gravitational field is 
assumed to be stationary, the chosen coordinate system is then a standard post-Newtonian gauge and the metric 
reduces to its usual form 



00 



1 



-W - 



21 



where W given by (|l]) reduces to 



0(6), {Go,} = ^^^^W^ + 0(5), G,, 



1 + ) 5., + 0(4) , (74) 



W{x) = U{x) + - 



p*{x') 



7+ 2 M 



(l-2/3)[/ + n- 



37— 
P 



d^x' . 



(75) 



and W is given by (48). As a consequence, for a clock moving with the coordinate velocity v, the quantity 1/m° is 
given by the formula 



,0 



1 



1 



W 



-V 



u" at \ 2 

from which it is easily deduced that 



7 + 



0(6) 



(76) 



u 







,0 



1 



1 



Wa-Wb 



-A 



^|(7 + 1){Wav\ - WbvI) + ]^{Wa - Wb) [Wa - T^b + «i - 4 + 2(1 



-P){Wa + Wb)] (77) 



1 



2 2 

■"aVb 



1 



v%-2{^ + 1){Wa-va-Wb-vb] 



0(6). 



It follows from ( [73| ) and (|7^) that the frequency shift Sv/v is given by 

^'A 



6v_ 

V 



- 1 = 



Vb 



(78) 



where {6v/v)c is the special-relativistic Doppler effect 



'5v 



1 



— = — Nab ■ {va- Vb) + ^ 



\'"A - - {Nab ■ {va - vb)) {Nab ■ vb) 
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{Nab ■ (va-vb)) 



\vl - ^vl + {Nab ■ vb) 



l^A - \'"Wb - \v% ~ {Nab ■ {va - vb)) {Nab ■ vb) Qw^ - + {Nab ■ vg) 
and {5v)/v)g contains all the contribution of the gravitational field, eventually mixed with kinetic terms 



(79) 
0(5) 



:{Wa-Wb) 



1 r 

73 



{Wa 



Wb) {Nab ■ {va ~ vb)) - I^a^ ■ va 



,(2) 



Vb 



+ + ^)iWAv\ - WbvI) + -{Wa- Wb) [Wa -Wb + 2(1 - I3){Wa + Wb) + vl 

,(2) 



(80) 



Va Vb 



-2 {Nab ■ {va - vb)) {Nab ■ vb)] + Nab ■ [(if ■ vb) {va - 2vb) + {l 
+ (Z^') - 2(7 + 1)Wa) ■ VA - {if - 2(7 + 1)Wb) ■ vb} + 0(5) . 

It must be emphasized that the formulas ( |76| ) and (^^ are valid within the PPN framework without adding special 
assumption, provided that P and 7 are the only non-vanishing post-Newtonian parameters. On the other hand, ( |80| ) 
is valid only for stationary gravitational fields. In the case of an axisymmetric rotating body, we shall obtain an 
approximate expression of the frequency shift by inserting the following developments in (pO|), yielded by (p3)-(p9|): 



;(2 



Im{xa,xb) +Ij2{xa,xb) 



' ^aVc^ = Is{xa,xb) H , 

• • > ^sVc^ = Is{xb,xa) H 

the function Is being now given by (|68|) written with ai — 0. Let us recall that the symbol 
contributions of the higher multipole moments which are neglected. 



Vg'/c = -Im{xb,xa) - Ij2{xb,xa) 



stands for the 



B. Application in the vicinity of the Earth 

In order to perform numerical estimates of the frequency shifts in the vicinity of the Earth, we suppose now that 
A is on board the International Space Station (ISS) orbiting at the altitude H = 400 km and that B is a terrestrial 
station. It will be the case for the ACES mission. We use = 6.37 x 10^ m and va — tb = 400 km. For the velocity 
of ISS, we take va = 7.7 x 10'^ m/s and for the terrestrial station, we have vb < 465 m/s. The other useful parameters 
concerning the Earth are: GM = 3.986 x 10^'' m^/s^, Ve = 6.378 x 10^ m, J2 = 1.083 x 10"^; for n > 3, the multipole 
moments J„ are in the order of 10"^. With these values, we get Wb/c^ ~ GM/c^tb = 6.95 x 10-^° and Wa/c^ ~ 
GM/c^TA =6.54x10^^". From these data, it is easy to deduce the following upper bounds: \Nab-va/ c\< 2.6x10"^ 
for the satelhte, | Nab ■ vb/c\< 1.6 x 10~^ for the ground station and | Nab ■ {va ~ vb)/c |< 2.76 x 10^^ for the 
first-order Doppler term. 

Our purpose is to obtain correct estimates of the effects in (^0|) with are greater than or equal to 10~^* for an 
axisymmetric model of the Earth. At this level of approximation, it is not sufhcient to take into account the J2- 
terms in {Wa — Wb)Ic^ ■ First, the higher-multipole moments J3, J4, ... yield contribution of order 10^^'^ in Wa/c^ ■ 
Second, owing to the irregularities in the distribution of masses, the expansion of the geopotential in a series of 
spherical harmonics is probably not convergent at the surface of the Earth. For these reasons, we do not expand 
{Wa - Wb)/^ in 

However, for the higher-order terms in (po|), we can apply the explicit formulas obtained in the previous section. 
Indeed, since the difference between the geoid and the reference ellipsoid is less than 100 m, Wb/(? may be written 
as II 

1 GM GMrlJ2, o , 1 . 

^^B = ^ + ^^(l-3cos^.) + -AI^., 

where the residual term AWb/c^ is such that | AW^/c^ |< 10^^^. At a level of experimental uncertainty about 10^^^, 
this inequality allows to retain only the contributions due to M, J2 and S in the terms of orders 1/c? and l/c"*. As a 
consequence, the formula (|o|) reduces to 

— =-{Wa-Wb) + -[ — ] +- — +••• 
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r 


\^ J M 






' s 



(81) 



where the different terms involved in the r.h.s. are separately explicited and discussed in what follows. 
Using the identity (r^ + r^)^ — R\b ~ 27'a^s(1 + '>T'A • ''t-b), it niay be seen that {Sv/v)^-^j is given by 



(3) 
M 



GMirA+rs) 



7 + 1 



TA - rs 



Nab ■ {va - vb) 



+ (7 + 1)- 



Rab ua-va + ub ■ Vb 



(82) 



rA + TB l + UA -riB 

The contribution of this third-order term is bounded by 5 x 10^^'* for 7 = 1, in accordance with a previous analysis 



C. Influence of the quadrupole moment at the order 

Defining the quantity Kab by 

Kab 



(ta - tbY 



tatb 

^(3) 



it is easily deduced from (^ and ( pO| ) that the term {5v/v) j^ in ( ^ ) is given by 

_ GM 

J-2 



2r, 



-J2 (Nab ■ {vA - Vb)) 



[3(fe • ua) 



1 



[3(fc-ns)' 



7 + 1 GMJ2rlirA+rB) 



1 



' A' B 



(Nab ■ (va-vb)) 



{l + UA- TIb)^ 

{k -UA + k- tlb) 



5 - Sua ■ nB + ^Kab 
l + UA -riB 



1 

Rab 



M(fc • ub)'^ + rsik ■ Ua)" 



TA + TB 

{ua -VA + nB- VB){k riA + k- ns) 



{i - Ua - ub + Kab) 



^7 - Ua - riB + 2Kab 



rA + TB 

u . „ 

{ua-va) [l-3(fc-n^)2] 



TA 

Rab 
tb 



I + ua-ub 
TA + rB{2 + nA- tib) 
ta + rB 



I \\^ 'ilv, .21 ^aC^ + riA ■ Ub) + TB 

[uB ■ Vb) [I - 3(fe • ns) J — 



+Rai 



/ UA ■ VA ^ riB ■ Vb 

\ TA TB 



l-(fc-ns)^ I 
-{Ua ■ Va) [ub ■ vb 



TA + TB 

[k ■ nA){k ■ ub) 

l-{k-nA)^ 



rB 



rA 



-2^{k-VA) 

TA 

r, Rab . 
-2 (fc • vb) 

rB 



, rA + 7-5(2 + UA ■ nB) , , 

k ■ riA [- k ■ Ub 

rA+TB 

rA{2 + nA ■nB) + rB 



k ■ UA + k ■ Ub 



rA + rB 



(83) 



One has | va/c |= 2.6 x 10^^, | vb/c \< 1.6 x 10"*^ and Kab 3.77 x 10"^. A crude estimate can be obtained by 
neglecting in ( p3| ) the terms involving the scalar products ub -vb and k-vs- Since the orbit of ISS is almost circular, 
the scalar product • va can also be neglected. On these assumptions, we find for 7 = 1 



1 



< 1.3 X 10" 



(84) 
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As a consequence, it will perhaps be necessary to take into account the 0(3) contributions of J2 in the ACES 
mission. This conclusion is to be compared with the order of magnitude given in Q without a detailed calculation. 
Of course, a better estimate might be found if the inclination i = 51.6 deg of the orbit with respect to the terrestrial 
equatorial plane and the latitude tt/2 — 9b of the ground station were taken into account. 



D. Frequency shifts of order 1/c* 



The term (Si^/i^)'-^^ in ( |8l| ) is given by 

(4) 



= (7 + 1) —vl 



M 



GM 2\ GM{rA~rB), 2 2 

1 '"A I 7. 



-{vl-vD 



1 / GM 



2 yrArs 
GM{rA + rs 



[{rA-rBf + 2{(3-\){r\~rl)\ 
2(7 + 1) ta- rs 



(85) 



tatb 
7 + 1 



l + UA-riB ta + tb 
Rab 



{Nab ■ (va-vb)) {Nab ■ vb) 



1 + ua ■ riB TA + tb 
,2 



{{uA ■ va) {Nab ■ vb) - {Nab ■ {va ~ 2vB)){nB ■ vb)) 



The dominant term (7 + l)GMvj^/rA in ( |85| ) induces a correction to the frequency shift which amounts to 10 . So, 
it will certainly be necessary to take this correction into account in experiments performed in the foreseeable future. 

The terms {5y /y)'^^^ is the contribution of the intrinsic angular momentum to the frequency shift. Substituting jsi] ) 
and (|6^ ) into (|80| ), it may be seen that 



^-^T = i.^s)A-{J's)B. 
/ S 



(86) 



where 



{^s)a = (7 + 1) — IH - \ — 

rW TB \\ + nA-nB 



"Tb 



k ■ {uA X Ub) 
{l + UA- TIbY 



VA + rB{2 + nA- tib) 



TA + TB 



TA + rB 



riA + riB 



■ k X UA 



(87) 



{^s)b = (7+1)^ 

^B 

_ fc ■ {uA X Ub) 
{l+UA- Ub)'- 



k X ua 



rA 



I + UA ■ riB TA + TB 
rA{2 + UA ■ TIb) + TB 

UA H ■ Ub 

TA + rB 



k X Ub 



In order to make easier the discussion, it is useful to introduce the angle between xa and xb and the angle ip 
between the plane of the photon path and the equatorial plane. These angles are defined by 

cos tp = Ua ■ nB , < tp < TT , k ■ {ua x ub) ^ sin ■0 cos ip , < ip < ir . 

With these definitions, it is easily seen that 

k ■ {ua X Ub) Ip 

— cos ir, tan — . 

l + riA- riB 2 



Let us apply our formulas to ISS. Due to the inequality vb/va < 6 x 10 , the term {J's)b '^^ ( pq ) may be neglected. 
From (|87|), it is easily deduced that 



\{^s)a\< (7 + 1)-^ 

^A 



rA\ 2 + 3 I tan -0/2 1 



rB 



1 1 + cos Ip I 



■ VA . 
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Assuming < ■)/' < 7r/2, we have (2 + 3 | tan -0/2 |)/ 1 1 + cos'0 |< 5. Inserting this inequaHty in the previous one and 
taking for the Earth GS/c'r\ = 3.15 x IQ-^^, we find 

< (7 + 1) X 10"^^ (89) 

Thus, we get an upper bound which is shghtly greater than the one estimated by retaining only the term hoiV^/c 
in {\n\}. However, our formula confirms that the intrinsic angular momentum of the Earth will not affect the ACES 
experiment. 



VI. CONCLUSION 



In this paper, we have shown that the world-function fl^XA, xb) constitutes a powerful tool for determining the 
coordinate time transfer and the frequency shift in a weak gravitational field. Our main results are established within 
the Nordtvedt-Will PPN formalism. We have found the general expression of fl{xA,XB) up to the order l/c^. This 
result yields the expression of the time transfer function T{tA,XA,XB) at the order 1/c"'. We point out that 7 and 
«! are the only post-Newtonian parameters involved in the expressions of the world-function and of the time transfer 
function. 

We have treated in detail the case of an isolated, axisymmetric rotating body, assuming that the gravitational field 
is stationary and that the body is moving with a constant velocity Vr relative to the universe rest frame. We have 
given a systematic procedure for calculating the terms due to the multipole moments in the world- function Q{xa, xb) 
and in the time transfer function T{xa, xb)- These terms are obtained by straightforward differentiations of a kernel 
function. We have explicitly derived the contributions due to the mass M , to the quadrupole moment J2 and to the 
intrinsic angular momentum 5' of the rotating body. 

Restricting our attention to the case where only /3 and 7 are different from zero, we have then determined the 
general expression of the frequency shift up to the order 1/c^. We have obtained the contributions of J2 at the order 
Our method would give as well the quadrupole contribution at the order 1/c^ in case of necessity. We have found 
the complete evaluation of the effect of the intrinsic angular momentum S, which is of order 1/c^. It is noteworthy 
that our formulas contain terms which have not been taking into account until now. 

Within the limits of our model, the formulas that we have established yield all the gravitational corrections to the 
frequency shifts up to 10~^* in the vicinity of the Earth. We have applied our results to the ACES mission. We have 
found that the influence of the quadrupole moment at the order 1/c^ is in the region of 10~^^. For the effect of the 
intrinsic angular momentum, we have obtained an upper bound which is greater than the currently accepted estimate 
but which remains three orders of magnitude less than the expected accuracy in an experiment like ACES. Finally, 
it must be noted that our results could be applied to the two-way time/frequency transfers. In particular, the 0(3) 
contributions of J2 to the two-way frequency transfers would probably deserve to be carefully calculated. 
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